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We use QCD sum rules to calculate the branching ratio for the production of the meson X (3872) 
in the decay B — > A(3872)7V assumed to be a mixture between charmonium and exotic molecular 
[cg][gc] states with J PC = We find that in a small range for the values of the mixing angle, 

5° < < 13°, we get the branching ratio B(B -> XK) = (1.00 ±0.68) x 10~ 5 , which is in agreement 
with the experimental upper limit. This result is compatible with the analysis of the mass and decay 
width of the mode J/4>(nn) and the radiative decay mode J/ipy performed in the same approach. 
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I. INTRODUCTION 

The A (3872) state has been first observed by the 
Belle Collaboration in the decay B + -> X(3872)K + -> 
J/ipir + Tr~ K + [l|, and was later confirmed by CDF, DO 
and BaBar This was the first state of an increas- 
ing number of candidates for exotic hadrons discovered 
recently. The current world average mass is mx = 
(3871.4 ± 0.6) MeV, and the width is T < 2.3 MeV at 
90% confidence level. BaBar Collaboration reported the 
radiative decay mode A(3872) — > yJ/ip [1, 3, which de- 
termines C = +. Further studies from Belle and CDF 
that combine angular information and kinematic proper- 
ties of the 7r + 7r _ pair, strongly favors the quantum num- 
bers J PC = 1++ or 2"+ % [|, Q. Although the new 
BaBar result favors the J PC = 2 ^ assignment es- 
tablished properties of the A (3872) are in conflict with 
this assignment . Therefore, in this work we will con- 
sider the A(3872) as being a J PC = 1++ state. BaBar 
Collaboration reported the upper limit of the branching 
ratio for the production in B meson decay [To| : 



B{B ± -> K ± X(3872)) < 3.2 x 10" 



(1) 



Recently, Belle Collaboration presented the most pre- 
cise measurement of the branching fraction B{B i 



A73872)fs: ± )S(X(3872)7J/V>) = (1.78±^±0.12)xl0- 6 
[llj . 

The decay modes of the AT (3872) into J/ip and other 
charmonium states indicate the existence of a cc in its 
content. However the attempts to classify the state in 
the charmonium spectrum have to deal with the fact that 
the mass of the X(3872) is not compatible with any of 
the possible candidates in the quark model [T||. An- 
other problem comes from the measurement of the de- 
cay rates of the processes A(3872) — > ,7 / 'rp ir + tt~ tt° and 
A(3872) — > J/ipn + Tr~, which are comparable Q, 



X -> J/V'TT+TT-TT 
X -^J/VJ-K+TT- 



= 1.0 ±0.4 ±0.3. 



(2) 
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that could indicate a strong isospin and G parity viola- 
tion, which is incompatible with a cc structure. 

The coincidence between the A (3872) mass and the 
D*°D° threshold: M(D*°D°) = (3871.81 ± 0.36) MeV 
[l3j . inspired the proposal that the A (3872) could be 
a molecular (D*°D° - D*°D°) bound state with small 
binding energy Il5j |. In particular, considering the 
A(3872) as an admixture of neutral and charged compo- 
nents of molecules, i.e. D°D*° and D + D*~, the strong 
isospin violation observed in Eq. (0) could be explained 
in a very natural way [Tfjl \l7\ . There is also a possibility 
that the observed ratio of the X decaying into J/tp + 2ir 
or 37t may not come from a large isospin breaking. In 
Ref. [l8| the isospin breaking is investigated in the dy- 
namical generation of the A as a molecular state, and it 
is found to be small. But, considering that the two pion 
and three pion states comes from the decays of p and 
oj mesons, the small isospin breaking is compensated by 
the larger phase space of the p meson, thus explaining 
the experimental data. 

Another interesting interpretation for the X (3872) is 
that it could be a compact tetraquark state (19h21| . 

The co-existence of both cc and multiquark compo- 
nents is subject of debate in many works, and it is sup- 
ported by some experimental data. In Ref. [22[ , a simula- 
tion for the production of a bound D°D*° state with bid- 
ing energy as small as 0.25 MeV, reported a production 
cross section that is an order of magnitude smaller than 
the cross section obtained from the CDF data. A similar 
result was obtained in Ref. [23[ in a more phenomeno- 
logical analysis. However, as pointed out in Ref. 24|, 
a consistent analysis of the D°D*° molecule production 
requires taking into account the effect of final state in- 
teractions of the D and D* mesons. 

Besides, the recent observation, reported by BaBar 
[H, of the decay A(3872) V ij{2S)j at a rate: 

B( B(X^7) 7) = 3 - 4 ± L4 ' [t is mUGh bi §g er than tllC 

molecular prediction r( ^^ ) ) 7) - 4 x 10~ 3 . 

In the framework of the QCD sum rules (QCDSR) the 
mass of the X were computed with good agreement with 
data, considering tetraquark [27} and molecular struc- 
tures [ItJ • The same success is not found in decay widths 
calculations. In Ref. [28| the decay width of the modes 
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J/ifj(mr) are calculated for the tetraquark structure, and 
the result is one order larger that the total width. In 
Ref. [1^, [3(| the QCDSR approach was used to study the 
X(3872) structure including the possibility of the mixing 
between two and four-quark states, where it was suc- 
cessfully applied to obtain the mass of the state and the 
decays widths for the modes J/ip(mr) and the radiative 
decay mode J/ipJ- This was implemented following the 
prescription suggested in [lH for the light sector. The 
mixing is done at the level of the currents and is ex- 
tended to the charm sector. In a different context (not 
in QCDSR), a similar mixing was suggested already some 
time ago by Suzuki [23j . Physically, this corresponds to 
a fluctuation of the cc state where a gluon is emitted 
and subsequently splits into a light quark-antiquark pair, 
which lives for some time and behaves like a moleculc-likc 
state 

The models for the quark structure can also be applied 
to study the production of the state in B decays. This 
subject is studied in different approaches in Refs. [32|-[34j]. 
In this work we will focus on production of the X(3872), 
usin g th e mixed two-quark and four-quark prescription of 
Ref. [29L [30l] to perform a QCDSR analysis of the process 
B ± X(3872)A' ± . 



II. THE DECAY B -> X(3872)A' 




The process B — > X(3872)A occurs via weak decay 
of the b quark, while the u quark is a spectator. The 
X meson as a mixed state of molecule and charmonium 
interacts via cc component of the weak current. In effec- 
tive theory, at the scale [i ~ m& <C mw, the weak decay 
is treated as a four-quark local interaction described by 
the effective Hamiltonian (see Fig. [T|): 



Gf 



C 2 (/x) 



02 



(3) 



the matrix element in two pieces: 

x {B(„)\J^\K( P '))(X(,,)\J""'>\0), (4) 

where p = p' + g, and the currents are 

J^ = ar M b, J<r> = cT>. (5) 

The matrix elements in Eq. are parametrized in the 
following way 



(X(q)\J^\0)=X w e;(q), 



(6) 



and 



{B(p)\J?\K(p')) = f + (q 2 )(p,+p^ + f-(q 2 )(p„ -*g . 

(7) 

The parameter Avk in ^ gives the coupling between the 
current and the X state. The form factors f±(q 2 ) 
describe the weak transition B — > K. Hence we can see 
that the factorization of the matrix element describes the 
decay as two separated sub-processes. 

The decay width for the process B ± -> X(3872)K ± is 
given by 

T(B -> XK) = ^^\V\ m l, m 2 K , m 2 x )\M\ 2 , (8) 

with X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. The 
invariant amplitude squared can be obtained from (U): 



\m\ 2 = ?f\VM 2 (c 2 + ^ 



x \{m 2 Bl m 2 K ,m 2 x )\ 2 w fl{q 2 )l 



(9) 



We will use QCD sum rules in order to determine the 
parameter \w and the form factor f+(q 2 ), and therefore 
we can obtain the width of the decay. 



III. TWO-POINT CORRELATOR 



The QCD sum rule approach [35M37| is based on the 
two point correlator: 

IV(g) - i [d 4 y e l ™(0\T{jX(y)J^}\0) , (10) 



where Vik are CKM matrix elements, C\(p) and C-x[)X) 
are short distance Wilson coefficients computed at the 
renormalization scale [i ~ 0(m\^). The four-quarks effec- 
tive operator is O2 = (cr At c)(sr M 6), with = 7^(1—75). 

The decay amplitude of the process is calculated from 
the Hamiltonian ([3]) , and it can be factorized by splitting 



where the current J^^ is defined in ([S]). For the X meson 
we will follow [3ll | and consider a mixed charmonium- 
molecular current as in Ref. [29L l30j . For the charmonium 
part we use the conventional axial current: 



fu 2) (x) = c Q (a;)7^75Ca(a;)- 



(11) 
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The D° D*° molecule is interpolated by |38H40j (with 
q = u): 



{la (z)75 Ca (x)c b (a) 7^5 {x) ) 



(g a (a;)7 AI c a (a;)cb(x)75gb(x)) 



(12) 



As in Rcf. [31j we define the normalized two-quark cur- 
rent as 



,•(29) 



1 



6^ 



(13) 



and from these two currents we build the following mixed 
charmonium-molccular current for the V(3872): 



J*{x) = MW 1 (x) + cos(fl)j W (x). (14) 

We will consider a small admixture of D + D*~ and 
£)-£)*+ components, then we have for the X current: 



Jrf (x) = cos a J' 1 (x) + shmJ^(x), 



(15) 



with J£(x) and J*{x) given by Eq.([T4"|). 

Considering the u and d quarks to be degenerate, i.e., 
m u = m e i and (uu) = (dd) , and inserting the currents (J5J) 
and (|15[) in the correlator we have in the OPE side of the 
sum rule 



n° PE (q) = (cos a + sin a) (sin 9 E$ (q) 
(«?) 



6^ 



coBfln^(g) 



(16) 



where 

n^ 4 (<z) = * | A e ^(0|T{J^(y)J KEc) (0)}|0) 
n^ 2 (g) = ijd 4 y e^y{0\T{J^(y)J u{5c) (0)}\0).(l7) 

(cc) 

The contribution from the vector part of the current Ji, 
vanishes after the integration is performed, hence these 
correlators are equaH except for a minus sign) to the ones 
calculated in Ref. |29| for the two-point correlator of the 
V(3872). 

On the phenomenological side the correlator is deter- 
mined inserting the intermediate state of the X: 



q 2 - m\ 



\jX\X( q ))(X{q)\jW\0) 



iX x X w I q^Qu 
r>2 _i_ ™2 [ 9nv 



X 



X 



(18) 



where q 2 = —Q 2 , and we have used the definition and 

(0\J?\X(q))=\ x e ll (q). (19) 

The parameter defining the coupling between the current 
J* and the X meson has been calculated in Ref. [29| , and 
its value is \ x = (3.6 ± 0.9) x 10~ 3 GcV 5 . 



In the QCDSR approach a Borcl transform to Q 2 — s- 
M 2 (Q 2 = —q 2 ) is performed to improve the matching 
between both sides of the sum rules. The Borel trans- 
form exponentially suppresses the contribution from ex- 
cited states in the phenomenological side of the sum rule. 
In the OPE side the Borel transform suppresses the con- 
tribution from higher dimension condensates [13] • After 
performing the Borcl transform we get in the structure 
9nu- 



X w X x e ^ = -(cos a + sin a) ( sin6»n 24 (M 2 ) 



6V2 



cos9W^(M 2 ) 



(20) 



This expression is analysed numerically to obtain the cou- 
pling parameter Xw- We perform the calculation using 
the same values for the masses and QCD condensates 
listed in [2j|, and in the same region in threshold pa- 
rameter so and Borcl mass M 2 that we have used in 
the mass and Xx analysis in Ref. (29|, ^/sq = 4.4 GcV, 
2.6 GcV 2 < M 2 < 3.0 GcV 2 . The mixing angles deter- 
mined in the same reference are: 



5° < 9 < 13° , a = 20° 



(21) 



Taking into account the variation in the Borel mass pa- 
rameter and the mixing angle 9, the result for the Xw 
parameter is: 



X w = (1.29 ±0.51) GeV 2 



IV. THREE-POINT CORRELATOR 



(22) 



The form factor of the B — s- K transition matrix ([7]) 
can be evaluated from the three point correlator: 

n M (p,j/) = J d\d'ye^'- x -P-y\0\T{J^(0)J K (x)jl(y)}\0) 

(23) 

In the OPE side we use the weak current Jj^ defined 
in ([5]) and the interpolating currents of the B and K 
pseudoscalar mesons: 

Jk{x) = iu a {x)~f 5 s a (x) , J B = i u a (x)j 5 b a (x) .(24) 

We work at leading order in a s and consider condensates 
up to dimension 5 and terms linear in the mass of the s 
quark. 

The phenomenological side of the sum rule is com- 
puted by inserting the intermediate states of the B and 
K mesons in the correlator: 



yrphen 



fBf K m 2 K m 2 B (/ + ft)(p + y) M + /_(t)g M ) 
m b (m s + m u ) (p 2 - m 2 B )(p' 2 - m 2 K ) 



(25) 
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with t 



(p ~ p') 2 ' using ([7]) and the following definitions: 



(0| J K \K(p')) = f K - 



K 



\J B \B(p)) = f B - 



(26) 



Performing a double Borel transform, P 2 — > M 2 and 
P' 2 — > M' 2 , and matching both sides of the sum rule, we 
get in the structure (p M + p'n) (with P 2 = —p 2 , P' 2 = 



p' 2 , Q 2 



-t): 



f B fKm 2 K m 2 B f+(t) 



-1 
4^2 



so rv-o m 2 

ds / dup pert (s, t, u)e~ a? 3 " - m 72 + -(ww)(m& + m s )e~~^ 

min 



mb(m u + m u ) An 2 J. . ./n ' 2 

(m 3 , + mjm s )(M 2 + M' 2 ) - M 2 (m b {M 2 + 2M' 2 + t)+ m s (M' 2 + 2M 2 + 1)) ) , 



rnZiutije m 



8M 4 M' 2 



where the perturbative contribution is given by 



p pert (s,t,«) 
x (2mbm 



u(2m b 



t + u) 



4A2 (s,t,u) 
s — u + f) + ((s — m 2 )(s — i + u) — 2su) 

2\ 



x (2?7ifcTO s — s + t — u) + (s + u — m6 )A(s, t, u) 



The integration limit s m j n is given by 



(28) 



(29) 



So = (m^+As) 2 and uq = (rriK+A u ) 2 are the continuum 
threshold parameters for the B and K respectively. Note 
that, after the double Borel transform, the contributions 
from the quark condensate and mixed condensates of the 
s quark are eliminated. 

We use the following relation between the Borel masses 
M 2 and M' 2 (H: 



M 



12 



0.64 GcV" 



M 



(30) 



A. Result for the form factor 

The sum rules are analysed numerically using the fol- 
lowing values for quark masses and QCD condensates, 
and for meson masses and decay constants [13, S3, 5I| : 

m b (m b ) = 4.7 GcV, m s = 0.140 GcV 
(qq) = -(0.23±0.03) 3 GeV 3 , 
m§ = 0.8 GeV 2 , 
m B = 5.279 GeV m K = 493.677 MeV 
f B = 0.170 GeV , f K = 0.160 GcV . (31) 



(27) 



I 

We use the value of the mixing angles a and 9 given in 
(f2"Tj) . For the continuum threshold parameters we take 
A s = A„ = 0.5 GeV. 

In Fig. [2] we show the plot for the form factor / + 
calculated in the sum rules from Eq. (|2T|) as a function 
of the Borel mass M 2 and the momentum transfer Q 2 . 
For the region M 2 > 20 GeV 2 shown in the plot the sum 
rule presents a good stability. 

In the following analysis, to determine the Q 2 depen- 
dence of the form factor, we choose the Borel mass within 
the region of stability around the mass of the B meson, 
26 GeV 2 < M 2 < 30 GeV 2 . We fit j27|) within the stable 
region by matching both sides of the sum rule. In Fig. 
|3] we show, through the dots, the QCDSR results for the 
form factor f+{Q 2 ) as a function of Q 2 . The numerical 
results can be well fitted by a monopolar parametrization 
(shown by the solid line in Fig. ([3])): 



U(Q 2 



(17.55 ±0.04) GeV 
(105.0 ± 1.76) GeV 2 + Q 1 



(32) 



For the decay width calculation, we need the value of 
the form factor at Q 2 = —m 2 Xl where mx is the mass of 



the off-shell X meson. Therefore we have: 



/+(Q 2 )lc3-- 



0.195 ±0.003. 



(33) 



B. The decay width B -> XK 

To determine the the decay width we insert, in the 
expression ([5]), the parameter \\y (f2"2"| obtained in the 
two-point sum rule calculation, and the value of the form 
factor /_)_, from the three-point sum rule, at the X meson 
pole (|3"3"|) . The branching ratio is therefore calculated 
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FIG. 2. The plot of the form factor /+ calculated as a function 
of Q 2 and M 2 from Eq. (27]). 
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FIG. 3. Momentum dependence of C(Q 2 ) for A s = A„ = 
0.5 GeV. The solid line gives the parametrization of the 
QCDSR results (dots) through l[3^l. 



dividing the result by the total width of the B meson 
Ttot = Ti/t b : 



B(B -> X('3872)K) = (1.00 ± 0.68) x 10 



-5 



(34) 



where we have used the mean life of the B meson tb = 
1.638 x 10" 12 s, the CKM parameters V cs = 1.023, V cb = 
40.6 x 10" 3 Hi, and the Wilson coefficients Ci_(/i) = 
1.082, C2 (u) = —0.185, computed at /x = m& and Ams = 
225 MeV j44[. The result (jM]) is in agreement with the 
experimental upper limit ([T]). 

For completeness we also compute the branching ratio 
for the X as pure cc and molecular states. We choose 
the mixing angle as 9 = 90° and 6 = 0° in Eq. (fi"4]) . 
and we get respectively for the pure molecule and pure 



B(B -> X mol K) = (0.38 ± 0.06) x 10" 



B(B -> X- CC K) = (2.68 ± 0.50) x 10" 



(35) 



(36) 



Comparing the results for the pure states with the one 
for the mixed state (|34[) . we can see that the branching 
ratio for the pure molecule is one order smaller, while the 
pure charmonium is larger. 

The result for the pure molecular state in Eq. (|35[) can 
be compared with the one from Ref. [33]. In this work 
the authors study the production of the X as a molecular 
state in the decay B — > XK through the coalescence of 
charm mesons. The calculation of the branching ratio 
is strongly dependent on the choice of parameters, and 
it is found to be of order 10~ 4 to 10~ 6 . In Ref. [Hj], 
the production of the X as a 2 3 Pi charmonium state 
is studied in pQCD, and their result for the branching 
ratio is 7.88^ 3 'y6 x 10 -4 , which is an order bigger than 
our result for the pure charmonium (|36[) . 



V. SUMMARY AND CONCLUSIONS 

We have presented a QCDSR analysis of the pro- 
duction of the X(3872) state considering a mixed 
charmonium-molecular current in the decay B — > XK. 
We find that the sum rules result in Eq. (|3~4"]) . obtained by 
using the factorization hypothesis, is smaller, but com- 
patible with the experimental upper limit. Since, it is 
known that non-factorizable contributions may play an 
important role in hadronic decays of B mesons [451 ] , our 
result can be interpreted as a lower limit for the branch- 
ing ratio. 

This result was obtained by considering the mixing an- 
gles in Eqs. ([T5j) and (|T4")) with the values a = 20° and 
5° < < 13°. The present result is also compatible 
with previous analysis of the mass of the X state and 
the decays into J/ipTT°Tr + Tr~ and J/ipTr + TT~ [2!|, and the 
radiative mode 'yJ/tp [3fj| . since the values of the mixing 
angles used in both calculations are the same. It is im- 
portant to mention that there is no new free parameter in 
the present analysis and, therefore, the result presented 
here strengthens the conclusion reached in Refs. 0, [3(| 
that the X(3872) is probably a mixture between a cc 
state and D°D*°, D°D*°, D+D*- and D~ D*+ molec- 
ular states. From Eq. (fT4"]) one may be tempted to say 
that the cc component of the state described by such cur- 
rent is dominant (~97%), like the conclusion presented 
in Ref. (29|. However, from a closer look at Eq. (|13l) . 
one can see that the cc component of our current is al- 
ready multiplied by a dimensional parameter, the quark 
condensate, in order to have the same dimension of the 
molecular part of the current. Therefore, it is not clear 
that only the angle in Eq. ([T4"]) determines the percentage 
of each component. To try to evaluate the importance of 
each part of the current it is better to analyse the results 
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obtained with each component, like the results presented 
in Eqs. and (|36D . From these results we see that 
the cc part of the state plays a very important role in the 
determination of the branching ratio. By the other hand, 
in the decay X — > and X — > J / '-07T+ ir~ ir° , the 

width obtained in our approach for a pure cc state is [2^] : 

T(X- CC -> J/i>rm) = , (37) 

and, therefore, the molecular part of the state is the only 
one that contributes to this decay, playing an essential 
role in the determination of this decay width. Also, for 
a pure cc state one gets: 

M XcE = (3.52 ± 0.05) GeV, (38) 



from where one sees again that the molecular part of the 
state plays a very important role in the determination of 
its mass. 

Therefore, although we cannot determine the percent- 
ages of the cc and the molecular components in the 
A (3872), we may say that both components are ex- 
tremely important, and that, in our approach, it is not 
possible to explain all the experimental data about the 
A (3872) with only one component. 
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